Fractal Stability Border in Plane Couette Flow 
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We study the dynamics of localised perturbations in plane Couette flow with periodic lateral 
boundary conditions. For small Reynolds number and small amplitude of the initial state the 
perturbation decays on a viscous time scale t oc Re. For Reynolds number larger than about 200, 
chaotic transients appear with life times longer than the viscous one. Depending on the type of the 
perturbation isolated initial conditions with infinite life time appear for Reynolds numbers larger 
than about 270-320. In this third regime, the life time as a function of Reynolds number and 
amplitude is fractal. These results suggest that in the transition region the turbulent dynamics is 
characterised by a chaotic repeller rather than an attractor. 
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For a variety of simple flows linear stability theory pre- 
dicts and experiments verify that the transition to tur- 
bulence proceeds by a sequence of instabilities at well de- 
fined critical values of the control parameter . Promi- 
nent examples include a layer of fluid heated from below, 
flow between rotating cylinders, liquid jets, and strati- 
fled fluids. The transition to turbulence in pipe flow |^] 
or plane Couette flow |^,|^ and to spiral turbulence in 
the flow between counter-rotating cylinders does not 
flt this pattern. If the Reynolds number is sufflciently 
large, turbulent dynamics can occur although the lami- 
nar profllc is still stable. A formal linear stability analysis 
predicts either no instability (as in plane Couette flow ) 
or one for Reynolds numbers larger than the ones where 
experiments begin to show turbulent behavior (e.g. plane 
Poiseuille flow Q). The nature of the transition is differ- 
ent from the cases with a linear instability. It depends 
on the size of the perturbation |^ , shows very strong in- 
termittency and has no sharply deflned stability 

border §. 

Such observations can be explained, if the bifurcations 
are sub critical, for then the new state extends to lower 
values of the Reynolds number and can be reached start- 
ing from flnite amplitude perturbations [p"2|~p^. Sim- 
ple models, based on interacting wave vector triads or 
truncations of Galerkin systems, are compatible with this 
view. However, it has not been possible to follow these 
upper branches all the way down to where they (dis- 
) appear and to identify a critical Reynolds number. 

A different approach has recently focused on the eigen- 
vectors of the linearised problem |]l^-|l8| . Very often the 
linearised hydrodynamic eigenvalue problem is not her- 
mitian, so that the eigenvectors are not orthogonal. As a 
consequence, even for (negative) eigenvalues, small per- 
turbations can be amplifled a lot, so that the neglect of 
nonlinear interactions no longer can be justified. This ex- 
plains how in the absence of linear instability the stability 
region around the laminar profile becomes small, perhaps 



irrelevantly small. One can imagine that because of this 
growth some of the turbulence seen in these systems is 
noise induced |Q. 

To shed some more light on the dynamics in these 
cases, we have studied the evolution of perturbations in 
a numerical model for plane Couette flow. The numer- 
ical method was designed specifically to allow accurate 
long integration times. The life time of these perturba- 
tions was mapped out in an amplitude-Reynolds number 
plane. This allowed to identify three dynamical regimes 
and to study the transitions between them. The results 
indicate that there is no sharp transition to turbulence: 
the landscape of life times has fractal features, with some 
isolated life times longer than the finite integration time. 
Thus, whether an initial condition leads to an evolution 
identified as turbulent depends in a sensitive way on ini- 
tial amplitude, Reynolds number and the experimental 
observation time. This is in qualitative agreement with 
experiments on plane Couette flow |^,|| and consistent 
with observations on pipe flow j^. In addition, the nu- 
merical results suggest that the turbulent state belongs 
to a transient repeller rather than a turbulent attractor 

The numerical model (see for more details) is based 
on an expansion of the velocity fleld in Fourier and Leg- 
endre polynominals. 



u(x, y, z,t)^Y. Uk,p(i)e'^'^"+'"'^ip(^) • 



(1) 



The advantage of Legendre polynominals compared to 
Chebyshev polynominals is that they give rise to evolu- 
tion equations which are energy conserving in the Eu- 
lerian, undriven limit. This is important for the long 
time evolution we want to study. For every velocity- 
component, we used a set of 37 Fourier modes on a hexag- 
onal grid and 16 Legendre-polynominals. The boundary 
conditions u{x, z, ±H/2) = and continuity equation are 
included within the Lagrange formalism of the first kind 
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PH . These restrictions finally leave 962 dynamically ac- 
tive degrees of freedom. 

The basic flow is normalised to have velocity ±Uo at 
z = ±H/2. The Reynolds number is defined as Re — 
UoH/2iy, with v the kinematic viscosity. Lengths are 
measured in units of H/2 and time in units of r = H/2Uq. 
The Navier-Stokcs equation for a perturbation u of the 
initial profile Uq UqzBx, 



u= -(u- V)Uo + (Uo- V)u 



-(u • V)u - -\/p + 

p Re 



Vu = 0, 



(2) 
(3) 



are replaced by nonlinear ordinary differential equations 
for the amplitudes Uk,p- It was verified that the time evo- 
lution could be followed sufficiently accurately for scaled 
times up to 10000. For h = 1 cm, the viscosity of water 
and a Reynolds number of Re = 400 this corresponds 
to about 9 min. In experiments a localised initial dis- 
turbance spreads and reaches the span wise walls within 
less than a minute |2^,|2^. The numerical calculations 
thus extend to a time region where in experiments the 
influence of the lateral walls are no longer negligible. 




FIG. 1. Sketch of the polodial vortex-ring {}^ used as initial 
condition. 



The initial state for our simulations was taken to be a 
poloidal vortex ring in the y-, z-plane, 



u = A curl curl (5 (a;, y)e 



(4) 



with a variable amplitude A. We approximate the S 
function by setting all real Fourier modes equal to one. 
As shown schematically in Fig. [l], the axis of the ring 
points in the e^j-direction. To remove excess energy for 
large modes we integrated this state for five time units 
at Re = 400 and used the resulting physically realistic 
initial condition in the subsequent studies. By this time 
the vortex ring has been rotated into the x-, y-plane and 
is similar to the one induced by the vertical water jets in 
the experiments of Berge et al. [||j2|,|2l . 



Three different flow regimes are distinguishable by 
their time series as shown in Fig. For low Reynolds 
number, the energy of the perturbation decays smoothly, 
with a bump related to the non-normality of the lin- 
earised equations of motion. The time when this maxi- 
mum is reached is in agreement with the estimate t„iax (x 
0.117 i?e in |2^. This is the first regime, dominated by 
the linearised equations. 
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FIG. 2. Evolution of the energy in the pertubation vs. time 
for three values of the Reynolds number. 



Around Re = 300 states show a qualitatively similar 
decay but only after an intermediate period of irregular 
dynamics. A life time can be defined by monitoring the 
z-component of the velocity until it decays below a cer- 
tain threshold. If that happens the perturbation lacks the 
driving and decays as described by the linearised equa- 
tions. In this second flow regime the life times for per- 
turbations are longer than estimated from the linearised 
equations. But one notes from Fig. ||that the life time of 
the disturbance is not a monotonic function of Re: The 
evolution at Re — 400 shows more fluctuations and vio- 
lent bursts but its life time is shorter than at Re ~ 300. 
In this third region, characterised by a non-monotonous 
variation of life times with amplitude and Reynolds num- 
ber, some initial conditions do not decay within the nu- 
merical integration time. The signal for the decaying and 
non-decaying states are dynamically similar in the fluc- 
tuating states, and there is no precursor that indicates 
the decay of the perturbation. 

The non-monotonous variations of life time with 
Reynolds number and amplitude occur for Re above 
about 350 and cover an increasing range in amplitude. 
On a global scale (Fig. ||), one notes a rather ragged 
landscape with a few peaks reaching up to the maximal 
life times followed in the numerical calculations. Because 
of the flnite observation times in experiments states are 
classified as decaying or turbulent according to the life 
times being above or below a certain threshold. Such a 
cut through the life time landscape shows isolated tur- 
bulent disturbances in the decaying region and isolated 
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decaying disturbances in the turbulent region. Such be- 
havior has been seen in experiments on constant mass 
flow through a pipe The non-monotonicity observed 
there may thus be of dynamical origin and not due to 
experimental limitations. 



the stable and unstable manifolds of the saddles so that 
trajectories trapped in this tangle need more and more 
time before they can escape and relax to the laminar 
profile. This is in line with previous observations in large 
systems and pipe flow that long lived transients 
can imitate a permanent turbulent state. Clearly a mea- 
sure of the thickness of the repeller or an estimate of the 
life time of transients as a function of Reynolds number 
would be highly desirable but is presently beyond our 
numerical or analytical abilities. 
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FIG. 3. Lifetime vs Reynolds number Re and energy of 
the initial state. The life time is normalised by Re, so that in 
the hnear viscous regime for small Re one expects a constant. 
Integration was stopped after a scaled time of 3000 units. 



The distribution of life times as a function of ampli- 
tude and Reynolds number has fractal features. Magni- 
fication of successive intervals in amplitude (Fig. ^) for 
fixed Reynolds number or of intervals in Re for fixed am- 
plitude reveals a self similar pattern of life times, without 
any obvious simplifications on smaller scales. The mag- 
nifications are reminiscent of life time pictures obtained 
in chaotic scattering p6|-^ . This analogy suggests that 
the long life times arise for initial perturbations which 
are close to the stable manifold of saddles in phase space. 
Then arbitrarily large life times should be possible and 
the 'turbulent' state would be supported by a repeller 
rather than an attractor. 

The magnifications in Fig. ^ give no indication that 
the peaks with large life times have a finite width, in 
accord with the repeller picture. The sensitivity of the 
life time to the initial state can be quantified in terms of 
the uncertainty dimension ||2^. For a given uncertainty 
in the amplitude e, we calculate the probability /(e) to 
find disturbance with life times greater than a given limit 
^max = 3000. For fractal sets, this probability decreases 
with decreasing e, with a scaling exponent a, /(e) oc e". 
The uncertainty dimension is then given by = 1 — a. 
For Re = 400 we find a scaling exponent a w 0.29 and 
an uncertainty dimension dg w 0.71. 

As the Reynolds number increases the number of long 
lived disturbances increases and they cover an ever in- 
creasing region in phase space and extend to smaller am- 
plitudes. Phase space becomes more densely filled with 
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FIG. 4. Successive magnifications of cuts through the 
life-time landscape of Fig. pi for Re — 400 and varying am- 
plitude (left column) and for fixed amplitude A — 100 and 
varying Reynolds number (right column). 



Further support for these observations is derived from 
a low dimensional system of ordinary differential equa- 
tions for a shear fiow where we can follow the dynamics 
more efficiently and can also study the nature of the re- 
peller in more detail In particular, we find that the 
repeller is built around a rapidly increasing number of 
stationary states which are born in saddle-node bifurca- 
tions above a critical Reynolds number. All states are 
unstable and the dimension of the repeller is rather high- 
dimensional. In the present flow the time intervals on 
the repeller are too short and the dimensions presum- 
ably to high to allow for a determination of its fractal 
dimension. An indication of its dimension is given by 
the number of unstable eigenvalues of the local lineariza- 
tion. In our simulation of plane Couette fiow at Re = 400 
this number lies between 30 and 80 positive eigenvalues. 
The maximum eigenvalue of the local linearization was 
about 0.35 in inverse time units. 

Candidates for repeller states in plane Couette fiow are 
the stationary states found by Nagata |0 and Clever and 
Bussc |3^j3^ . These states appear for Reynolds numbers 
above about 125. For larger Reynolds numbers we have 
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found further states, in accord with the observations on 
the model. But there remains a discrepancy between the 
Reynolds number for the occurrence of stationary states 
(about 125) and the one where the dynamics shows tur- 
bulence. However, the latter depends strongly on the 
initial conditions. For the ones presented above, the long 
lived states occur for Re > 300, but if one takes poloidal 
rings with instead of they occur for Re > 270 al- 
ready. This strong dependence on the nature of the initial 
states indicates that the stable manifold of the repeller 
are curled up in the high-dimensional phase space. As the 
Reynolds number increases so does the region explored 
until it eventually overlaps with the initial state. 

To summarise, we have presented evidence for a fractal 
transition region to long lived turbulent states in plane 
Couette flow. These turbulent states are transient and 
belong to a repeller. The density of long life times in- 
creases with Reynolds number. This scenario might be 
relevant for other shear flows without linear instability 
and perhaps for some systems with a subcritical transi- 
tion. 
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